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Abstract. The interaction of hydrogen with many transition metal sur- 
\^ ' faces is characterized by a coexistence of activated with non-activated 

, paths to adsorption with a broad distribution of barrier heights. By per- 

0^ ' forming six-dimensional quantum dynamical and classical molecular dy- 

namics calculations using the same potential energy surface derived from 
^ ^1 ab initio calculations for the system H2/Pd(100) we show that these fea- 

5^ ' tures of the potential energy surface lead to strong steering effects in 

^0 , the dissociative adsorption dynamics. The adsorption dynamics shows 

j^-^ ■ only a small isotope effect which is purely due to the quantum nature of 

I hydrogen. 

^ ' 1 Introduction 

. It is a long-term goal in surface science to understand catalytic reactions oc- 

I curing at surfaces. Obviously, the single steps of these often rather complicated 

. processes are more effectively studied at simple systems. The dissociative ad- 

sorption of molecules on surfaces is one of the fundamental reaction steps oc- 
. curing in catalysis. This establishes its technological relevance and importance. 

In particular, the dissociative adsorption and associative desorption of hydro- 
gen on metal surfaces has served as a benchmark system, both experimentally 
and theoretically (see, e.g., Refs. ||, ||, Q and references therein). Since the 
mass mismatch between hydrogen and a metal substrate is rather large, the 
crucial process in the dissociative adsorption for these particular systems is the 
conversion of translational and internal energy of the hydrogen molecule into 
translational and vibrational energy of the adsorbed hydrogen atoms. If in ad- 
dition no surface rearrangement occurs upon adsorption, the substrate degrees 
' of freedom can be neglected and the dissociation dynamics can be described in 

terms of potential energy surfaces (PES) which take only the molecular degrees 
of freedom into account. 

The PES for the dissociative adsorption of a diatomic molecule neglecting the 
substrate degrees of freedom is still six-dimensional. These PESs now become 
available by elaborate density- functional calculations H, ||, 0, |[ ||- However, in 
order to understand the reaction dynamics one has to perform dynamical cal- 
culations on these potentials. Because of its light mass hydrogen has to be de- 
scribed quantum mechanically. Only recently it has become possible to perform 
dynamical calculations of the dissociative adsorption and associative desorption 
where all degrees of freedom of the hydrogen molecule are treated quantum 
mechanically . These calculations have established the importance of high- 
dimensional effects in the reaction dynamics. 



X 



For example, molecular beam experiments of the dissociative adsorption of 
H2 on various transition metal surfaces like Pd(lOO) [0, Pd(lll) and Pd(llO) 
W(lll) W(IOO) H, W(100)-c(2x2)Cu |l| and Pt(lOO) Q 

revealed that the sticking probability in these systems initially decreases with 
increasing kinetic energy of the beam. High-dimensional quantum dynamical cal- 
culations have shown that steering effects can cause such an initial decrease in 
the sticking probability |l^; it is not necessarily due to a precursor mecha- 
nism, as was widely believed. 

So far we have studied, besides the steering effect the influence of the 
molecular rotation and orientation jl^, ^ and vibration on the sticking 
probability of H2/Pd(100). In this paper we present a comparison of quantum 
and classical dynamics for the dissociation of hydrogen on Pd(lOO) and inves- 
tigate isotope effects. In the next section the theoretical background will be 
introduced before the results of the dynamical calculations will be discussed. 



2 Theoretical background 

The potential energy surface of Il2/Pd(100) has been determined using density- 
functional theory together with the generalized gradient approximation (GGA) 



21 and the full-potential linear augmented plane wave method p2| , p3| . Ab initio 
total energies have been evaluated for more than 250 configurations | M] an d have 
been parametrized in a suitable form for the dynamical calculations ||10| . 

Figure 1 shows a cut through the PES of Il2/Pd(100), where the most 
favourable path towards dissociative adsorption is marked by the dashed line. For 
this path there is no energy barrier hindering dissociation, i.e., the adsorption is 
non-activated. However, the majority of pathways towards dissociative adsorp- 
tion has in fact energy barriers with a rather broad distribution of heights and 
positions, as the detailed total-energy calculations showed [Q, i.e. the PES is 
strongly anisotropic and corrugated. This has important consequences, as will 
be shown below. 

The quantum dynamics is determined in a coupled-channel scheme within the 
concept of the local reflection matrix (LORE) |24|. This numerically very stable 
method is closely related to the logarithmic derivative of the solution matrix and 
thus avoids exponentially increasing outgoing waves which could cause numerical 
instabilities. The reported calculations, which take all degrees of freedom of the 
hydrogen molecule into account, are still only possible if all symmetries of the 
scattering problem are utilized. 

The classical trajectory calculations are performed on exactly the same PES 
as the quantum dynamical calculations. The equations of motions are numeri- 
cally integrated with the Bulirsch-Stoer method with a variable time-step [^. 
The energy conservation is typically fulfilled to 0.1 meV. The sticking proba- 
bility is determined by averaging over a sufficient number of trajectories. The 
exact number of trajectories depends on the specific initial conditions and ranges 
between 1,815 and 18,330. 




Fig. 1. Contour plot of the PES along 
a two-dimensional cut through the 
six-dimensional coordinate space of 
H2/Pd(100). The inset shows the ori- 
entation of the molecular axis and the 
lateral H2 center-of-mass coordinates. 
The coordinates in the figure are the H2 
center-of-mass distance from the sur- 
face Z and the H-H interatomic dis- 
tance d. The dashed line is the opti- 
mum reaction path. Energies are in eV 
per H2 molecule. The contour spacing 
is 0.1 eV. 
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As far as the CPU time requirement is concerned, it is a wide-spread believe 
that classical methods are much less time-consuming than quantum ones. This 
is certainly true if one compares the computational cost of one trajectory to a 
quantum calculations. However, if one is interested in averaged quantities like 
sticking probablities, then in classical molecular dynamics calculations one has 
to average over many trajectories corresponding to different initial conditions. 
Quantum mechanics does this averaging automatically. A plane wave describing 
the incident beam hits the surface everywhere, and a j = rotational state con- 
tains all molecular orientations. Thus for the results presented here the quantum 
method is even more time-efficient than the classical calculations, in particular 
if one considers the fact, that in a coupled-channel method the sticking and 
scattering probabilites of all open channels are determined simultaneously. 



3 Results and Discussion 

Figure 2 presents six-dimensional quantum dynamical calculations of the 
sticking probability as a function of the kinetic energy of a H2 beam under nor- 
mal incidence on a Pd(lOO) surface and five-dimensional calculations for D2. In 
addition, the results of a H2 molecular beam experiment are shown . Quan- 
tum mechanically determined sticking probabilities for hydrogen at surfaces with 
an attractive well exhibit an oscillatory structure as a function of the incident 
energy ^ reflecting the opening of new scattering channels and 
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Fig. 2. Sticking probability versus kinetic energy for a liydrogen beam under nor- 
mal incidence on a Pd(lOO) surface. Experiment (H2): circles (from ref. the- 
ory: six-dimensional results for H2 molecules initially in the rotational and vibra- 
tional ground state (dashed line) and with an initial rotational and energy distri- 
bution adequate for molecular beam experiments (solid line), and vibrationally adia- 
batic five-dimensional results for D2 molecules initially in the rotational ground state 
(dash-dotted line). 



resonances [p6[ |27| . These structures are known for a long time in He and H2 
scattering [p8[ and also in LEED For H2/Pd(100), however, measuring these 
oscillations is a very demanding task. They are very sensitive to surface imper- 
fections like adatoms or steps and therefore hard to detect ||3^ . Since we do 
not focus on these oscillations here, for the solid line in Fig. 2 we have assumed 
a velocity spread of the incoming beam typical for the experiment pH so that 
the oscillations are smoothed out. 

The initial decrease of the sticking probability found in the experiment is well- 
produced in the quantum dynamical calculations. The high sticking probability 
at low kinetic energies is caused by a steering effect: Slow molecules can very 
efficiently be steered to non-activated pathways towards dissociative adsorption 
by the attractive forces of the potential. This mechanism becomes less effective 
at higher kinetic energies where the molecules are too fast to be focused into 
favourable configurations towards dissociative adsorption. This causes the initial 
decrease of the sticking probability. If the kinetic energy is further increased, the 
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Fig. 3. Probability for dissociative adsorption versus kinetic translational energy for a 
H2 beam under normal incidence on a clean Pd(lOO) surface for non-rotating molecules. 
The solid line corresponds to the six-dimensional quantum dynamical results assum- 
ing an energy spread typical for beam experiments. The molecular dynamics results 
(H2: circles, D2: diamonds) have been obtained purely classical, i.e. without including 
zero-point energies in the initial conditions. The inset shows an enlargement of the 
results at low energies. 



molecules will eventually have enough energy to directly traverse the barrier 
region leading to the final rise in the sticking probability. 

In Fig. 3 we compare the averaged quantum mechanical sticking probabihty 
for H2 with the results of classical trajectory calculations for H2 and D2. The in- 
set shows an enlargement of the results at low energies. The molecular dynamics 
calculations have been performed without any zero-point energies in the initial 
conditions. This allows us to truely differentiate between classical and quan- 
tum effects in the dynamics. First of all, the classical results do not show any 
oscillatory structure revealing that the oscillations are purely due to quantum 
mechanics. Furthermore, at low energies the classical results fall almost exactly 
upon the averaged quantum results. This shows that steering is a general concept 
and not restricted to quantum or classical dynamics. We attribute the difference 
between the classical and quantum results at the medium energy range to the 
influence of zero-point effects in the multi-dimensional interaction potential 132] . 



At very high energies, where zero-point efects play no role, quantum and classical 
results again agree. 

Since the averaged quantum and the classical sticking probability agree in 
the low-energy regime, we can use classical trajectories to describe the steering 
mechanism. This is done in Fig. 4, where snapshots of two typical trajectories 
are shown. The initial conditions are chosen in such a way that the trajectories 
are restricted to the xz-plane. 

The left trajectory illustrates the steering effect The initial kinetic 

energy is Ei = 0.01 eV. Initially the molecular axis is almost perpendicular to 
the surface. In such a configuration the molecule cannot dissociate at the surface. 
But the molecule is so slow that the attractive forces can reorient the molecule 
so that it can follow a non-activated path towards dissociative adsorption. 

In the case of the right trajectory, the initial conditions are the same as in 
the left one, except that the molecule has a higher kinetic energy of 0.12 eV. 
Due to the anisotropy of the PES the molecule also starts to rotate to a con- 
figuration parallel to the surface. However, now the molecule is so fast that it 
hits the repulsive wall of the potential before it is in a favorable configuration to 
dissociative adsorption. It is then scattered back into the gas-phase rotationally 
excited. 




Fig. 4. Snapshots of classical trajectories of hydrogen molecules impinging on a Pd(lOO) 
surface. The initial conditions are chosen in such a way that the trajectories are re- 
stricted to the xz-p\ane. Left trajectory: initial kinetic energy Ei = 0.01 eV. Right 
trajectory: same initial conditions as in the left trajectory except that the molecule 
has a higher kinetic energy of 0.12 eV. 



Finally we like to discuss the isotope effects in the dissociation of hydrogen 
on Pd(lOO). Fig. 3 shows that in clasical dynamics there is practically no isotope 
effect between H2 and D2 in the sticking probability. As far as the low-energy 
regime is concerned, this seems surprising at a first glance, since D2 is more inert 
than H2 due to its higher mass. However, one has to keep in mind that at the 
same kinetic energy D2 is slower than H2, so that there is more time for the 
steering forces to redirect the D2 molecule. This compensating effect has also 
been found theoretically for the dissociation of hydrogen at W(IOO) |ll8| . 

In Fig. 2 we have also plotted the sticking probability of D2 according to five- 
dimensional quantum dynamical calculations. Due to its higher mass the energy 
spacing between the quantum levels is smaller for D2 than for H2 . Therefore much 
more eigenfunctions in the expansion of the wavefunction have to be taken into 
account in the coupled-channel calculations for D2 than for H2. This makes a six- 
dimensional quantum treatment of D2 not feasible at the moment. However, we 
have already shown that the results of five-dimensional vibrationally adiabatic 
quantum calculations, where the molecules are not allowed to make vibrational 
transitions, are very close to the full six-dimensional results for the dissociation 
of hydrogen on Pd(lOO) j2^. Hence it is reasonable to compare five-dimensional 
results for D2 with six-dimensional results for H2. The quantum dynamical stick- 
ing probabilities of D2 are slightly smaller than those of H2 . Since no such isotope 
effect is observed in the classical calculations (Fig. 3), this small difference has to 
be a quantum mechanical effect. We attribute it to the larger vibrational zero- 
point energy of H2 which can effectively be used to traverse the barrier region 
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4 Conclusions 

In conclusion, we reported a six-dimensional quantum and classical dynamical 
study of dissociative adsorption of H2/Pd(100). We have shown that the initial 
decrease of the sticking probability is due to a steering mechanism which is 
operative both in quantum and classical dynamics. The adsorption dynamics 
shows only a small isotope effect which is purely due to the quantum nature of 
hydrogen. Our results establish the importance of a high-dimensional dynamical 
treatment in order to understand reactions at surfaces. 
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